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1. Introduction. 

The possibility that higher- derivative corrections should be added to the Einstein-Hilbert (EH) la- 
grangian of general relativity in order to obtain a better behaved theory has been often considered. Among 
the various possibilities, a preminent role is played by the so-called Gauss-Bonnet (GB) terms. 

GB lagrangians were introduced in [1] as the only possible generalization of the EH lagrangian in higher 
dimensions that gives rise to field equations which are second order in the metric, linear in the second 
derivatives, and divergence free. Another important property is that GB corrections do not introduce any 
new propagating degrees of freedom in the spectrum of gravity [2]. However, since they vanish in lower 
dimensions unless nonminimally coupled to scalar fields, they are mainly useful in the context of higher- 
dimensional gravity, and especially Kaluza-Klein theories [3-5]. It must however be mentioned that GB 
contributions also appear in the low-energy limit of string theories, [2,6], and may also play an important 
role in the context of the braneworld scenario [7] . 

The introduction of GB terms in the action of Kaluza-Klein theories allows spontaneous compactification 
of higher dimensional models without the need of introducing external fields. For example, GB models admit 
ground states in the form of the direct product of two maximally symmetric spaces [4]. They also have 
interesting applications in higher-dimensional cosmology [5]. 

In order to explore further the physical implications of the dimensional reduction of higher dimensional 
models of gravity including GB corrections, it is interesting to study the existence of black hole solutions 
of the dimensionally reduced theory, with compact internal space. In the case of pure Einstein gravity, this 
investigation was performed in [8] , where it was shown that the only solution of physical interest is the four- 
dimensional Schwarzschild metric with fiat internal space. In the GB case, one may expect the existence of 
a greater variety of solutions, and in particular also black holes with anti-de Sitter asymptotics. 

Some black hole solutions of the Einstein-GB field equations are already known in difl[erent physical 
situations, as spherical symmetry in higher dimensions [9] or GB-scalar coupling in four dimensions [10,11]. 
In these cases it results a modification with respect to the Einstein case of the short-distance behaviour of 
the solutions near the singularities, but also asymptotic or global properties of the black hole may be altered. 

In this paper, our aim is to classify all solutions of the Einstein-GB system taking the form of a direct 
product of a four- dimensional spherically symmetric black hole with a maximally symmetric internal space. 
Since a general discussion would be too involved, we shall limit ourselves to the case of six dimensions, 
where the only relevant GB correction is quadratic in the curvature and has the form <S = TZ^'^^'^TZnv pa — 

A powerful technique for investigating this topic is the study of the phase space of the solutions of 
the field equations. This method has been used for example in the Einstein case [8]. As mentioned above, 
when a GB term is added to the action, the field equations are still second order, and linear in the second 
derivatives, but no longer quadratic in the first derivatives. This fact gives rise to several technical problems. 
In particular, the potential of the dynamical system is no longer polynomial, but presents poles for some 
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values of the variables [11]. 

The result of our investigation is that physically relevant black hole solutions exist that are not asymp- 
totically flat. 

Let us start by considering the (n + 4)-diniensional action 

1 = j ^ S^+^^x (7e("+^) + a<S("+^)), (1.1) 

where TZ^"-~^^'> is the curvature scalar and iS^"^''-' the quadratic GB term of the manifold and a is a coupling 
constant of dimension [L]'^. 

We want to perform a dimensional reduction which casts the metric in the form of a direct product of 
a four-dimensional manifold with an n-dimensional space of constant ciirvature. whose size is parametrized 
by a scalar field 0. In general, contrary to the Einstein case, it is not possible to find an ansatz for the 
metric of the Einstein-GB system that completely disentangles the scalar field d from the curvature in the 
dimensionally reduced action, except when the internal space is flat. Therefore we maintain the usual ansatz 

dsln+4) = e-"*cis^4) + e'^gi;^dx'^dx\ (1.2) 
where rfs^_^^ is the line element of the four-dimensional spacetime and g'^^ is the metric of the n-dimensional 
maximally symmetric internal space, with TZ^^ = Xig^^ ■ The action is dimensionally reduced to 



j d^x [(1 + 2aXi e-2^)7e(^) + a e"'^^^^) + Ana e"'^a^t^ V-^V".; 
nin + 2) ,„2_. _1o^^v.-20^ t^.,2_ nin + 2){n^ + n--i) 



I 

+ l^llil^^ - - 2n - 12)aAi e'^-^J (V<^)^ - ' ^"'^^ ' a e"^{Vct>y 

+ Xi e-("+2)'^ + („ _ 2)(n - 3)aA2e-("+^)^l . (1.3) 

The dimensionally reduced action contains the Einstein and the GB terms nonminimally coupled to 
a scalar field, and a standard kinetic term and a potential for the scalar field. In addition, one has a 
non-standard quartic correction to the kinetic term and a coupling between the Einstein tensor Q^i, and 
derivatives of the scalar field. Of course, both these terms yield second order field equations. The action 
(1.3) displays some similarity with the string effective action studied in [11], but contains additional terms. 

In the following discussion, it is important to fix the possible ground states for the model. These are 
taken to be the direct product of a four dimensional and an n-dimensional maximally symmetric space, i.e. 
TZ'ttvpa = ^eigj^pgil^ - g^p}g\^p ), Tlp'Jpa = ^iiglTjgi'^ - glTJgl'p')- Substituting this ansatz into the field 
equations derived from (1.1), one obtains 

a[(n - l)(n - 2)(n - 3)(n - 4)Aj + 2AKl + 2n(n - l)KeM + {n - l){n - 2)Ai + 12Ae = 0, 

a[n(n - l)(n - 2)(n - 3)A'^ + 12n(n - l)AiAe] -t- n(n - l)Ai + 6Ae = 0. (1.4) 

The system always admits the solution Ag = Aj = 0, as in the Einstein case, but one can also obtain solutions 
with uoiu auisliiug (•ur\-atur(\ uamely de Sitter or anti-de Sitter.* Consequently, black hole solutions of (1.1) 



* We are only interested in black hole with asymptotic regions, so we shall not consider the de Sitter case 
further. 



may have anti-de Sitter behavior at spatial infinity. In the following we shall concentrate on the case n = 2. 
Eq. (1.4) then admits a solution Ag = Aj = i.e. AdS x H'^ for a > 0, ot dS x S'^ for a < 0. 



2. The dynamical system. 

Let us consider the case n = 2. For the four- dimensional metric we adopt the spherically symmetric 
ansatz [8] 

dsfi^ = -e^'^dt'^ + a-^e^^-^-'df + e"^-^" Qijdx^ dx\ (2.1) 

where v, Q and a as well as cf) are functions of ^ and gab is the metric of a two-dimensional maximally 
symmetric space, with TZij = XeQij- Of course, in the case of physical interest, Ag > 0. 
It is then convenient to define new variables 

X = 2C - 1/ - <A, r? = 2C - - 2</. (2.2) 

Substituting the ansatz (1.2), (2.1) into the action, performing some integrations by parts, and factoring out 
the internal space, the action can be cast in the form 



1= -Stt j d'^xi^a 



6x'' + 3C'2 + 3?7'2 - 8x'C' - 8x'r/' + K'i] - -(A^e^C + A.e^") 



^{ri - x')(4C' + 3ry' - + a(C' - x')(3C' + V - Sx^AiS^" 

- o\c; - x'W - X')(llx" + + + 7Cv' - 13x'C' - 13x'r?') - ^eA^ |. (2.3) 

As usual, the action (2.3) does not contain derivatives of a, which acts therefore as a Lagrangian 
multiplier enforcing the Hamiltonian constraint. Another relevant property of (2.3) is that, in spite of the 
presence of the higher derivative GB term, it contains only first derivatives of the fields, although up to the 
fourth power, and therefore gives rise to second order field equations. A further interesting property is that, 
due to our choice of variables, the action is invariant under the interchange of ( and r/. 

One can now vary (2.3) and then write the field equations in first order form in terms of the new 
variables, 

W = x', X = C, Y = rj', U = e^, Z= e^, V= e", (2.4) 

which satisfy 

U' = WU, Z' = XZ, V' = YV. (2.5) 
Varying with respect to a and then choosing the gauge (7 = 1, one obtains the Hamiltonian constraint 

E=P'^+ XeZ^ + XiV"^ 

+ % [AeAiZV^ + Ae^2(y - W)A + \iV'^{X - W)B - 3(X - W)i^ - W)^^] = 0, (2.6) 



where 



P2 = ^ _|_ _ _ _^ ^2 ^ -^^^2 _^ 4^2 _^ 4y2 _^ _ 13^^ - UWY, 

A = 4X + 3Y -5W, B = SX + AY- 5W. 

Variation with respect to x, ^ and gives rise to the other field equations 

{2a 
jj^ [XeZ'^{2X + AY- 5W) + XiV'^{4X + 2Y- 5W) + 22W^ 

- 2X^ - 2Y^ - 36W'^X - 36W^Y - UX'^Y - 12Y^X + 17WX^ + VJWY"^ + A4XYW] |' 

= ^ [ - XeXiZ^V^ + XeZ^{Y - W)A + XiV^{X - W)B - {X - W){Y - W)C^] , (2.7) 

X' + 2Y' - 2W' + [AeZ2(2X + 2Y- 3W) - {X - W){WW'^ + 2X'^ + 5Y'^ + 6XY 

- 9WX - UWY - XiV^)] I' = XeZ^ + ^ [X^V^{X - W)B - {X - W){Y - W)C'^] , (2.8) 

{4rv 
jj^ [XiV'^{2X + 2Y- 2.W) -{Y - W){IQW'^ + SX^ + 2Y'^ + &XY 

- UWX - 9WY - XeZ'^)] }' = XiV'^ + ^ [XeZ'^{Y - W)A - {X - W){Y - W)C'^] . (2.9) 

In the variables (2.4), the problem takes the form of a six-dimensional dynamical system, subject to a 
constraint. Notice that the function E defined in (2.6) is a constant of the motion of the system (2.5), 
(2.7)-(2.9), whose value vanishes by virtue of the Hamiltonian constraint. Since the system is obviously 
symmetric for V —V, Z —Z, U —U, we shall only consider positive values of these variables. 

The Einstein limit 

In the Einstein limit, a = 0, one recovers the results of [8]. We summarize them in terms of the variables 
introduced above: when a = 0, the dynamical system reduces to eqs. (2.5) and 

2X' + 2Y' -3W' = 0, X' + 2Y' -2W' = XeZ"^, 2X' + Y' - 2W' = XiV'^ , (2.10) 

subject to the constraint 

E = p2 + Ae^^ + AiF^ = 0. (2.11) 

The physical trajectories lie on the four-dimensional hyperplane E = 0. Moreover, the system is independent 
of the variable U, and one may restrict the analysis to U = 0. It is evident that 2X + 2Y — 3W is a constant 
of the motion for the system (2.10) and therefore one of the variables, say W, could be eliminated, but we 
keep it for comparison with the GB case. 

The critical points at finite distance correspond to the small radius limit of the solutions. They lie on 
the surface Zq = Vo = Pq = 0, but only points with Xq = Yo = Wo correspond to regular horizons, while the 
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others give rise to naked singularities. The eigenvalues of the linearized equations around the critical points 
are 0(3), Xo,Yo,Wq. 

Since we are interested in solutions with asymptotic regions, we are led to study the phase space at 
infinity, which corresponds to the large radius limit of the solutions. This can be investigated defining new 
variables 

I X Y U Z V 

*=w' ^=w' y=w^ "=w^' '=w^ ''=w- 

In terms of these variables, the field equations at infinity are then obtained for i ^ 0, and read 

i = -2{v'^ + z^)t, X ^ + 2v^ - 2{v^ + z^)x y = 2z'^ + - 2{v^ + z^)y, 

u = {1 - 2v^ - 2z'^)u, z={x-2v^ -2z'^)z, v = {y - 2v^ - 2z^)v, (2.13) 

where a dot denotes t d/d(_. 

The critical points at infinity are found at = and 

a) AjWg = XsZq = 0, X ^ xq, y ^ j/g, with 3a:;o + Sy^ + 4a;gyg - 8xo - 8yo + 6 = 0. 

b) X,vl = 0, XeZ^ = 1/4, xa = 1/2, yg = 1- 

c) AeZg = 0, Xivl = 1/4, xo = 1, = 1/2. 

d) X,v^ = Xez^ = 3/16, xo = yo = 3/4. 

Points a) are the endpoints of the hypersurface V = Z = 0, points b) of the hypersurface V = 0, points c) 
of the hypersurface Z = 0. Exact solutions with endpoints b), c) and d) are discussed in the appendix. 

The eigenvalues of the linearized equations around the critical points and their degeneracy are: 

a) 0(5), 1, Xo, yo- 

b,c) -1,-1(5), i (2). 

rl^l -3 _3/9\ 1 3±vT5 
^) 2' 4 V^/' 4' 8 ■ 

The asymptotic behavior of the solutions can be deduced from the location of the critical points at 
infinity [8] . Excluding points a) that do not correspond to physical trajectories, one has, in terms of a radial 
variable r: 

b) ds'^ dt'^ + dr'^ +r'^dn\, e^"^ ~ const. 

c) ds'^ r^di^ + r^dr^ + r^dO^, e^"^ ~ r^. 

d) ~ -r dt^ + dr"^ + r'^dO?^, e'^'^' ~ r. 

Here, we have denoted with dfl\ the metric of a unitary 2-sphere, and with dOg that of a flat 2-plane. The 
solutions ending at points b) are asymptotically flat, while the others have more exotic behavior. 

The phase space portrait is the following: solutions with regular horizons start at = Vg = 0, 
Xq = Yo = Wo = a, for some value of the parameter a, and end at points b) if = 0, or d) if Aj > 0. These 
last solutions, however, decompactify for r — > oo, since e^*^ diverges in such limit. Also cylindrical solutions 
with Ae = exist, which end at points c). 



3. The Gauss-Bonnet phase space 

As discussed in section 1, in the GB case eqs. (1.4) admit the ground state solution Ag = —l/2a, 
Ki = — 3/lOa in addition to flat space, and therefore black holes with anti-de Sitter asymptotic behavior 
may be expected if a > 0. The phase space of the system can be studied by the same methods used in the 
Einstein case. As usual in the presence of Gauss-Bonnet terms, some care must be taken because of the 
poles sX U = Q. The limit [/ must be therefore taken, when necessary, at the end of the calculations. 

Equations (2.7)-(2.9) must be solved for the variables A', Y' and W in order to put the system in its 
canonical form. One can then find the critical points at finite distance by requiring the vanishing of the 
derivatives of the fields. As in the Einstein case, they lie on the hypersurface Uq — Za = Va — 0. However, 
in the GB system, the other variables must satisfy the constraint Wq = Aq = Iq, or Xq = ^^^Wo, 
Yq = ^^^Wq. Only the first case corresponds to regular horizons. In that case the eigenvalues of the 
linearized equations near the critical points are identical to those found in the Einstein limit. 

The critical points at infinity are obtained by writing the dynamical systems in terms of the variables 
(2.12) and requiring the vanishing of their derivatives as t — > 0. We find the following points: 







2/0 








a) 


1 


1 











b) 


1/2 


1 





1/4 





c) 


1 


1/2 








1/4 


e) 


2/3 


1 


2/9 





-1/15 


f) 


1 


2/3 


2/9 


-1/15 





g) 


1 


1 


2 


-1 


-1 


h) 


4/5 


4/5 


6/25 








i) 


2/3 


1 





1/3 





1) 


1 


2/3 








1/3 


m±) 


4±v'5 
R 


4=FV'5 












In case a), the curve of the previous section reduces to a single point. Furthermore, the points d) have 
disappeared, except in the limit a — > (i.e. t oc). The points with uq = are attained by taking the 
limit uo ^ at the end of the calculation. It is also evident that points e)-h) exist only if a > 0. 

It is interesting to notice that the location of the critical points at infinity is very similar to that of the 
pure Einstein system with a cosmological constant, investigated in [12], except for the presence of the new 
points i), 1) and m±). It seems therefore that one of the the main ingredients in fixing the structure of the 
phase space at infinity, and hence the asymptotic behavior of the solutions, is the relative dimension of the 
terms in the action. For a detailed discussion see [13] . 

Prom the eigenvalues and the eigenvectors of the linearized equations, one can deduce the nature of the 
trajectories attracted by the various critical points at infinity. The eigenvalues of the linearized equations 
near the critical points and their degeneracy are listed in the following table, together with the nature of the 
trajectories attracted, for W > Q. 
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Eigenvalues (with degeneracy) 


Trajectories attracted 


a) 


0(5) 1(5) 




b) 


-1(3) -1 ^(2) 


> A,- = 




-I (5), -1 4 (2) 


A^ = A, > 


c) 


W ^\ 1 1±a/i1/3 

-l(^), -2, -i, 


any Ae, A^ < 


f) 




Ae < 0, any A^ 


g) 


1, 2(^), 1, '^Y'^' 


Ae < 0, Aj < 


h) 


-1(5), -2,-i(^) 


any A,., A, 


i) 


,2' 3' 3 v^; 


Ae > 0, A, = 


1) 


~3' ~3' 3 ('^) 


Ae = 0, Ai < 


m+) 


2^0-1 n 1 2±3y5 
3 V^/' 3 ' 15 


Ae > 0, Ai = 


m_) 


2 co\ p 1 2±3\/5 

3 V'^/' 3 ' in 


Ae = 0, Ai < 



The points a) do not attract any trajectory from finite distance. 

The critical points b)-c) generalize those found in the Einstein case, and have the same asymptotic 
behavior. For what concerns the other points, 



e) 


ds^ ~ 


-r'^di^ + r-^dr^ + r^dO^, 


e2<^ 


^ const. 


f) 


ds^ ~ 


-r'^de + dr'^ + r'^dn'^_, 


e2^ 


rvj 7-2 _ 


g) 


ds^ ^ 


-r'^dt'^ + r-'^dr'^ + dfl'i, 


e2<^ 


^ const. 


h) 


ds^ ^ 


-r'^dt^ + r-^dr'^ + r^dfig. 


e2^ 


~ r. 


i) 


dS^r^ 


-r''de + dr'' + r''dn\, 


e2^ 


const. 


1) 


ds''^ 


-r^dt"^ +r'^dr'^ +r'^dfll, 


e20 




m±) 


ds^ ^ 




e2^ 





Here dfl'^ denotes the metric of a two-dimensional space H"^ with constant negative curvature. 

Of particular interest are the solutions that end at the critical point e), which arise for positive a. These 
asymptote to the exact ground state solution AdS^ x H'^, cited previously, that in the present coordinates 
takes the form 

ds^ = -(— + l]df+(— + l] \r'+r^dQl, 6^* = — . (3.1) 
\2a J \2a J 3 

Also interesting is the solution g), that asymptotes the exact solution AdS^ x H'^ x H^. Its four- dimensional 
section is analogous to a Bertotti-Robinson metric. The other solutions have less common behavior. 

The phase space structure can be summarized as follows: solutions with regular horizon start at the 
points U = V = Z = 0, X = Y = W and terminate at one of the critical points listed above, depending on 
the values of a, Ae and A,. 

Prom the Kaluza-Klein point of view, the relevant solutions are those with Ae > and e^'^ asymptotically 
constant. These are the asymptotically flat (Schwarzschild-like) solutions with flat internal space b), the 
asymptotically anti-de Sitter solutions with negatively curved internal space e), and possibly the more exotic 
the solutions i) and m+), again with flat internal space, whose nature is however not very clear. 
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4. Conclusions 

Higher- dimensional models of gravity naturally admit Gauss-Bonnet terms in the lagrangian. Our study 
has shown that the compactification of the simplest model admits black hole solutions displaying a variety 
of asymptotic behaviors. In particular, besides asymptotically flat and asymptotically anti-de Sitter black 
holes, also physically reasonable solutions with asymptotic behavior of type i) are present. 

It turns out that the phase space of the model is quite similar to that of pure Einstein gravity with 
a cosmological term [12]. It would be interesting therefore to consider the effect of adding a cosmological 
constant to our model. This topic is currently under study [13]. 
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Appendix A 

In the pure Einstein case, exact solutions corresponding to the vanishing of Xi or Ag where obtained in 
[8] for generic spacetime dimensions. In our six-dimensional setting, the solutions take the following form: 
for Aj = one obtains of course the Schwarzschild metric with constant scalar field, 

ds^ = -(1 - 2M/r)dt'^ + {1 - 2M/r)-^dr'^ + r^dn^, e^'^ = const. 

For Ae = 0, one has instead a solution of the form 

ds^ = -r2(l - 2M/r)df +r^il- 2M/r)-^dr^ + r^d^l, e^^ = . 

In our case, a family of exact solutions with asymptotic behavior d) can also be obtained, if one makes 
the ansatz r? = C- In fact, in this case, the field equations reduce to 

4C" - 3x" = 0, 3C" - 2x" = e^^, 

subject to the constraint 6%'^ + 10^''^ — 16x'C' + 2e^^ = 0. Integrating the first equation, one obtains 
x' = 4(C' — c)/3, and hence = Ae^^^"'^^'/^ for constant A and c. Substituting in the second equation, one 
obtains 

C" = 3e2«, 

which is solved by 

^ _ 2ae< 

Regular black hole solutions satisfying the constraint are obtained for c = a/4. For r] = the metric 
functions of (2.1) are related to our variables by 
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Using these relations with A = l, defining r = J e^^d^, ro = 2a/3, and substituting in (2.1), one finally 
obtains 

j-i/d r — ro 

or, in different coordinates, 

= (l - ^) + i (l - ^) " + = 
A more familiar expression can be obtained by writing the metric in its six-dimensional form: 

ds' = - (i - dt' + (i - 5 ) " df^ + Y idn^ + dnl). 

This is a variant of the well known six-dimensional Tangherlini metric, where the 4-sphere is replaced by the 
direct product 5^ x 5^. 

Appendix B 

It may be interesting to write down some special exact solutions of the Einstcin-GB system corresponding 
to the possible asymptotic behavior associated with the different critical points at infinity. The properties 
of these solutions are more transparent in their six-dimensional form, in the Schwarzschild-like gauge 

ds^ = -e^^dt" + e-2^dr2 + e^^dO^ ^ ^2a^^2_ 

In this gauge is evident the presence of a symmetry for the interchange of p and a, that follows from the 
specific compactification considered. This entails a duality between points b), c) and e), f). 

b) ds"^ = -de + dr'^ + r'^dn\ + dnl. 

c) ds"^ = -de + dr'^ + dnl + r'^dn\. 

e) ds'' = -{^ + l)dt'+(^ + iy\r^+r^dn\ + ^dn''_, or 

= -^de + ^dr''+ r^dQl + ifi dnl. 

f) = _ ri ^ 1^ ^^2 + (l^: + l) dr^ + l^dQl+ r^dn%, or 

ds^ = -^ de + f« dr^ + ifi dnl + r'^dnl 

g) = - - m) dV" + (i^ - m)" dr^ + 2adn'i + 2adn'i . 

h) ds'' = -^ de + ^ dr^ + r^dfi^ + ^2^Q2_ 
The parameter m is an arbitrary constant. 

We were unable to find exact solutions in the other cases. 
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